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Abstract: Using a criterion proposed by Salpeter and standard solid-state 
physics, we have determined conditions for the occurrence of the plasma-solid 
transition. Possible astrophysical applications are discussed. 

Introduction 

Phase transitions of various kinds have been, and still are, actively con- 
tributing to important changes occurring in the world around us. The aim 
of the present paper is a contribution to the study of the conditions for the 
occurrence of the so-called "plasma-solid (PS)" transition. 

Simple physical reasoning shows that the PS transition must certainly 
occur in various astrophysical settings, such as in proto-planetary and proto- 
stellar clouds. On the one hand, it is widely known that the universe contains 
different kinds of plasmas. On the other hand, it is equally widely known that 
solid objects also exist in the universe. This implies that there must exist a 
physical region where a transition between the two regimes takes place. 

In our preliminary work on the PS transition [1], we have determined the 
conditions for its occurrence in two simple idealized systems: a pure Fermi- 
Dirac and a pure Bose-Einstein gas. The object of the present calculation is 
again a determination of the conditions for the occurrence of the PS tran- 
sition, but this time on a more realistic footing. Our starting point is the 
criterion for the occurrence of the PS transition (although it was not called 
by that name) proposed by Salpeter [2]. 

Salpeter's paper was devoted to a thorough discussion of a zero-tempera- 
ture plasma. In its third part, he considered a system of positive ions of 
given charge and mass, rigidly fixed in the nodes of a perfect crystal lattice, 
and went on to estimate the zero point energy of the ions. He showed that 
the behavior of such a system can be described by the ratio 
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where E ZiP denotes the zero point energy of the ions and Ec is the Coulomb 
energy. According to the analysis presented in [2], a PS transition occurs for 
/ = 1. Calculations reported in [2] were estimates, and strictly valid only 
for T = K. Our calculations reported in the following section go beyond 
those assumptions in at least two aspects: (i) we have not used estimates 
but exact calculations, and (ii) our results take into account the influence of 
temperature. 



Calculations 



Finding a general expression for the energy of a real solid that takes into 
account most (if not all) of its characteristics is a formidable problem in solid 
state physics (for example [3] , [4] ). Various approximations to the complete 
solution exist; one of them is the Debye model of a solid. Within the Debye 
model, the energy per mole of a solid is given by the following expression [3]: 
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where N denotes the number of elementary cells in a mole of the material, 
n the number of particles per elementary cell, u(v) — jf — the static lattice 
energy per cell, 6 the Debye temperature, ks the Boltzmann constant, and 
in the following we will use the convention fc B = 1- The second integral in 
Eq. (2) can be solved as [5] 
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under the condition n >z 1, \x\ -< 2tt. Here, the symbol B 2 k denotes 
Bernoulli's numbers. It follows from Eq. (2) that the energy per particle of 
a solid within the Debye model is given by 

E P = 9T^I (4) 
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Using Eq. (3) in Eq.(5), one finally gets 
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which implies that 
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or, expanding explicitly 



E p = 3T 
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The energy of a Fermi-Dirac gas of particles of mass m and spin s enclosed 
in a volume V is given by [6] 
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All the symbols in this equation have their usual meaning, in particular, 
g — 2s + 1, and F 3 / 2 is a particular case of a Fermi integral 
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Eq. (10) can be developed as [7] 
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where T and q are the gamma function and Riemann zeta function, respec- 
tively. 
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Expanding to fourth order and inserting the temperature dependence of 
the chemical potential, one gets the following result for the energy per particle 
of the electron gas 
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The chemical potential of the electron gas at T = K is denoted here by 
= An 2 / 3 , with A = (3vr 2 ) 2/3 h 2 /2m. 

Introducing B = (2 1 / 2 /5)gm 3 / 2 / (tt 2 ^ 3 ) simplifies the notation in Eq. (12). 



Discussion 

Salpeter's ratio defined in Eq. (1) can be formed using Eqs. (8) and (12). 
The result obtained this way can be expressed as follows 
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Imposing the condition / = 1 on Eq.(13), it becomes possible to determine 
the values of the parameters of a PS transition. Solving Eq. (13) for the Debye 
temperature 9, and introducing in the solution the density dependence of the 
chemical potential at T = K one obtains the following expression 
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This equation links the relevant parameters of a plasma (number density and 
temperature) with those of the solid (Debye temperature) that can condense 
from it. We have thus obtained the equation of state of a system undergoing 
a PS transition. Apart from being interesting from the point of view of pure 
statistical physics, this result can find astrophysical applications in studies 
of dust and gas clouds. For recent examples of observational studies of such 
clouds, see for example [8] [9]. 

The number density and the temperature of an astrophysical cloud can 
be determined from observations using the methods of plasma physics, such 
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as the analysis of spectral line broadening, or in some cases, in-situ measure- 
ments of interplanetary gas and dust. Once these values are known for a 
given cloud, the Debye temperature of the solid material that can condense 
from it is given by Eq. (14). The Debye temperature is a unique character- 
istics of every solid, which gives hope that the method discussed here could 
lead to a possibility to determine the chemical composition of a condensing 
proto-planetary system. Such a determination would of course be only indi- 
rect, because after calculating the value of 9 from Eq. (14), one would have 
to rely on an identification of the material (or a mixture of materials) based 
on the corresponding calculated value. 

Consider Eq.(14) in the limiting case T — * 0, /io — * 0. Astronomically, 
this limit corresponds to an interstellar cloud sufficiently far away from any 
bright star. Developing Eq. (14) into series in the chemical potential and 
retaining just the first two terms, one gets the following expression for the 
Debye temperature 

9 = KxT^Vo V2 + K 2 T^^ 2 (15) 

The numerical values of the constants in Eq. (15) can be determined by 
developing Eq. (14) into a series. By doing this, it follows that 
^ 1 =(2/3) 1 / 2 5 1 / 6 7 1 / 4 7T 2 / 3 J B 1 /6 an d K^I^I^-^-^ihB)- 1 ^. 

Here, the constant B is the one defined after Eq. (12). For sufficiently small 
values of the chemical potential, Eq. (15) can be simplified into 

9 = KiT^/io 1/2 (16) 

Inserting the values of all the constants which occur in Eq. (16), one finally 
arrives at 

9 = 50.8T 7/6 n~ 1/3 (17) 

This expression has two possible applications. It can be used for the calcula- 
tion of 9 if the temperature and the number density of the cloud are known. 
Typical values of the temperature and density of the molecular clouds in our 
Galaxy are 20 < T[K] < 60 and 10 3 < n(cm' 3 ) < 10 5 . Inserting in Eq. (17) 
n = 2000cm -3 and T = 43 K gives 9 = 320 K, which is nearly the exper- 
imentally known value of the Debye temperature of the chemical element 
magnesium. This result could be interpreted as indication that in a molecu- 
lar cloud, grains of Mg can condense under these densities and temperatures. 
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In a similar way, the pair of values n = 10 4 cm~ 3 , T = 120 K would lead to 
6 = 628 K. Such a value is only slightly higher than the experimental value 
for the element Si. Note that this last result has a direct significance for the 
interpretation of observations, because emission from silicate grains has in- 
deed been detected [10]. Turning the argument around, if solid particles are 
observed in a cloud, and if their chemical composition can be determined, it 
becomes possible to calculate the Debye temperature from the principles of 
solid state physics. If, in addition, temperature is known from spectroscopy, 
one can determine the value of the chemical potential, and finally, the num- 
ber density of a cloud from Eq. (16). Further work along these lines is in 
progress and will be discussed elsewhere. 
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